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ABSTRACT 

An accurate method to measure the abundance of high-redshift galaxies consists in the 
observation of absorbers along the line of sight toward a background quasar. Here, we 
present abundance measurements of 13 z > 3 sub-Damped Lyman-a Systems (quasar 
absorbers with H i column density 19 < log N(HI) < 20.3 cm~^) based on the high 
resolution observations with VLT UVES spectrograph. These observations more than 
double the metallicity information for sub-DLAs previously available at z > 3. This 
new data, combined with other sub-DLA measurements from the literature, confirm 
the stronger metallicity redshift evolution than for the classical Damped Lyman-a 
absorbers. Besides, these observations are used to compute for the first time the frac- 
tion of gas ionised from photo-ionisation modelling in a sample of sub-DLAs. Based on 
these results, we calculate that sub-DLAs contribute no more than 6% of the expected 
amount of metals at 2,^ 2.5. We therefore conclude that even if sub-DLAs are found 
to be more metal-rich than classical DLAs, they are insufficient to close the so-called 
"missing metals problem" . 

Key words: galaxies: abundance - galaxies: high-redshift ~ quasars: absorption lines 
- quasars: PSS J0118-I-0320, PSS J0121-I-0347, SDSS J0124-I-0044, PSS J0133-t-0400, 
BRI J0137-4224, BR J2215-1611, BR J2216-6714. 



1 INTRODUCTION 

Damped Lyman-a systems (hereafter DLAs) seen in absorp- 
tion in the spectra of background quasars are selected over 
all redshifts independent of their intrinsic luminosity. They 
have hydrogen column densities, log A'^(HI) > 20.3. DLAs 
are also contributors to the neutral gas, rini, in the Uni- 
verse at high redshifts and it is from this gas reservoir that 
the stars visible today formed (Wolfe et al. 1995). Further- 
more, DLAs offer a direct and accurate probe of elemental 
abundances over > 90% of the age of the Universe. 

Recently, much attention has been drawn to the sub- 
Damped Lyman-Q Systems. These systems have H I col- 
umn density 19 < log N(HI) < 20.3 cm~^ and were first 
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coined by Peroux et al. (2003a). In a series of paper, our 
group has analysed a unique homogeneous sample of sub- 
DLAs all observed at the same resolution with the same 
instrument on the Very Large Telescope (VLT). The first 
two papers are based on the ESO/UVES archives (Paper 
I: Dessauges-Zavadsky et al. 2003; Paper II; Peroux et al. 
2003b) while more recent studies are from our own obser- 
vational programmes (Paper III: Peroux et al. 2005; Paper 
IV: this paper) leading to a new sample of high-redshift sub- 
DLAs observed under the same conditions. These works as 
well as other recent studies are suggesting that high metal- 
licities can be found more easily in sub-DLAs than in clas- 
sic DLAs, especially a low-redshift (e.g., Pettini et al. 2000; 
Jenkins et al. 2005; Peroux et al. 2006a; Prochaska et al. 
2006; Peroux et al. 2006b). In the present paper, we study 
a new sample of high-redshift sub-DLAs to better constrain 
their metallicities at z > 3. In addition, our data are used 
in combination with others from the literature to compute 
the fraction of ionised gas in a sample of sub-DLAs and re- 
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Table 1. The sample of 13 high-redshift sub-DLAs for which 
abundance studies have been undertaken. 



Quasar 






^abs 


log A^Hi 


PSS J0118- 


f0320 


4.230 


4.128 


20.02±0.15 


PSS J0121- 


f0347 


4.127 


2.976 


19.53±0.10 


SDSS J0124+0044 


3.840 


2.988 


19.18±0.10 








3.078 


20.21±0.10 


PSS J0133- 


f0400 


4.154 


3.139 


19.01±0.10 








3.995 


19.94±0.15 








3.999 


19.16±0.15 








4.021 


19.09±0.15 


BRI J0137- 


-4224 


3.970 


3.101 


19.81±0.10 








3.665 


19.11±0.10 


BR J2215- 


1611 


3.990 


3.656 


19.01±0.15 








3.662 


20.05±0.15 


BR J2216- 


6714 


4.469 


3.368 


19.80±0.10 



liably compute the contribution of sub-DLAs to the global 
metallicity. 

The paper is structure as follows. In the second sec- 
tion, we present the methodology and results of the deter- 
mination of the chemical content of these high-redshift sub- 
DLAs. In the third section, the total abundances including 
results from detailed photo-ionisation models are given and 
the redshift evolution of the metallicity is presented. Finally, 
the last section describes the contribution of sub-DLAs to 
the global metallicity in the context of the so-called missing 
metals problem. 



2 A SAMPLE OF 13 Z > 3 SUB-DLAS 

Table 1 summarises the absorption redshift and A'hi column 
densities for each of the 13 z > 3 sub-DLAs which constitute 
the sample under study. The data acquisition and reduction 
are described in Peroux et al. (2005). In the case of PSS 
J0133-I-0400, our own data have been supplemented by a 
3000 sec exposure with setting 860 (ESO 74.A-0306, PL: 
Valentina D'Odorico) and a 5200 sec exposure with setting 
540 (ESO 73.A-0071, PL: Cedric Ledoux). 

In this section, details on each of the studied sub-DLAs 
are provided. The metals associated with the Lyman lines 
were searched over the spectral coverage available. More 
than the 40 transitions most frequently detected in high 
column density quasar absorbers have been systematically 
looked for. The metal column densities are then determined 
by fitting Voigt profiles to the absorption lines. The fits were 
performed using the minimisation routine fitlyman in 
MIDAS (Fontana & Ballester 1995). 

(i) PSS J0118-I-0320 (zom=4.230, Zabs=4.128, 
logiVm=20.02±0.15): 

This system is at the high end of the sub-DLA definition 
with log A'^Hi=20.02±0.15. Many metal lines are detected at 
Zabs=4.128, some of which are clearly saturated. The fit is 
done using simultaneously Fe ii A 1608, Fe ii A 1144, Si ii 
A 1526, Si II A 1304, S ii A 1259 and S ii A 1253 to derive 
the redshifts and 6 parameters of the 5 components. The fit 
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Figure 1. Fit to the low-ionisation transitions of the 2abs=4.128 
logAfHi=20.02±0.15 sub-DLA towards PSS J0118+0320 (see Ta- 
ble 2). In this and the following figures, the zero velocity corre- 
sponds to the absorption redshift listed in Table 1, the vertical 
dashed lines correspond to the fitted components and the men- 
tion 'Ly forest' means that the metal line falls in the Lyman-a 
forest of the quasar spectrum while 'Tellurics' indicates a region 
contaminated by telluric lines. 



is shown in Figure 1 and the matching parameters are pre- 
sented in Table 2. In this and the following tables velocities 
and b axe in km s"'^ and N's are in cm~^. Lower limits for 
C II and O I are derived from the saturated lines. The Zn ll 
A 2026 line is covered by our spectrum but is situated in 
a region too noisy to allow to derive any meaningful upper 
limit. 

Both Si IV and C iv high-ionisation doublets are detected. 
The redshifts and b of the 9 components are derived from 
a simultaneous fit of C IV A 1548, C iv A 1550 and Si iv A 
1393 (Si IV A 1402 is affected by telluric hues). The resulting 
fit is overplotted on Si iv A 1402 for consistency check. The 
fit is shown in Figure 2 and the matching parameters are 
presented in Table 3. 

(ii) PSS J0121+0347 (zcm=4.127, 2abs=2.976, 
logiVm = 19.53±0.10): 

The low-ionisation transitions in this system are well fit- 
ted with one component, the redshift and b of which are 
fixed by a simultaneous fit of Si ll A 1260, O I A 1302, C ll 
A 1334, Fe ll A 1608 and Al ll A 1670. The first four lines 
are situated in the Lyman-a forest. This explains the likely 
blending which affects them. Therefore, measurements from 
Si III, O I and C ll are considered as upper limits. The fit 
is shown in Figure 3 and the matching parameters are pre- 
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Table 2. Parameters fit to the low-ionisation transitions of the Zabs=4.128 log 7Vhi=20.02±0.15 sub-DLA towards PSS J0118+0320. 



z 


b 


log N(Fe II) 


log N(Si II) 


log N{C II) 


log N(0 i) 


log N(S II) 


4.127836 


4.61±0.14 


13.21±0.11 


13.60±0.02 


>15.35 


>13.95 


13.41±0.11 


4.128319 


30.06±0.21 


13.57±0.12 


14.14±0.02 


>14.58 


>15.15 


13.94±0.11 


4.128678 


3.11±0.14 


13.71±0.11 


14.69±0.04 


>17.28 


>16.37 


13.75±0.11 


4.129264 


3.30±0.18 


13.43±0.12 


13.42±0.02 


>15.64 


>14.02 


13.02±0.13 


4.129627 


6.87±0.25 


13.12±0.11 


12.73±0.01 


>12.60 


>12.00 


11.89±0.12 



Table 3. Parameters fit 
of the Zabs=4.128 logArHi = 
J0118+0320. 



to the high-ionisation transitions 
:20.02±0.15 sub-DLA towards PSS 



z 


b 


log N(C iv) 


log N(Si IV) 


4.126680 


29.70±0.26 


13.49±0.12 


13.17±0.12 


4.127345 


9.90±0.13 


13.19±0.12 


12.47±0.12 


4.128024 


22.20±0.31 


13.30±0.11 


12.71±0.22 


4.128656 


11.10±0.12 


13.30±0.33 


12.56±0.12 


4.129083 


21.60±0.23 


13.59±0.11 


13.04±0.11 


4.129787 


3.20±0.02 


13.04±0.12 


12.30±0.12 


4.130161 


8.70±0.04 


13.50±0.11 


12.97±0.12 


4.130446 


2.60±0.01 


13.08±0.12 


12.65±0.12 


4.130711 


9.00±0.08 


13.31±0.11 


12.84±0.12 







~ Feii_i eoe 










Ly for&3-t~~^ 






Cll_1 334- 


Ly forest 






OI_ 1 302 


Ly forest 






-Alll_1 GVO , 





1.0 L 



















CIV_1 550 
































CIV_1 54S 
































SilV_1 402 










A 















































\ 








SilV_1393 












A 





Figure 2. Fit to the high-ionisation transitions of the Zabs=4.128 
log AfHi= 20.02±0.15 sub-DLA towards PSS J0118+0320 (see Ta- 
ble 3). 



sented in Table 4. AI ill A 1854 and Al ill A 1862 are covered 
but not detected. An upper limit is derived for the Al abun- 
dance: log N(A1 iii)<11.77 at 4o-. 

Interestingly, in this system the high-ionisation transitions 
is partially fitted with the same velocity component than for 
the low-ionisation transitions. The C iv and Si iv doublets 
are fitted simultaneously. C iv is heavily blended and there- 
fore leads to an upper limit, while Si iv is nicely fitted from 
the Si IV A 1393 line. The fit is shown in Figure 4 and the 
matching parameters are presented in Table 5. 



Figure 3. Fit to the low-ionisation transitions of the 2abs=2.976 
logAfHi=19.53±0.10 sub-DLA towards PSS J0121+0347 (see Ta- 
ble 4). 



Table 5. Parameters fit to the high-ionisation transitions 
of the 2abs=2.976 log AfHi=19.53±0.10 sub-DLA towards PSS 
J0121-I-0347. 



z 


b 


log N(C iv) 


log N(Si iv) 


2.975965 


11.90±0.08 


<13.50 


13.10±0.10 


2.976579 


16.90±0.08 


<13.47 


12.95±0.10 



(iii) SDSS J0124-f0044 (2cm=3.840, 2abs=2.988, 
logiVm = 19.18±0.10): 

This system has low-ionisation lines which cover a rather 
large velocity range and the profiles are well characterised by 
two separate clumps at ~ —100 km s~^ and ~ -1-30 km s~^. 
A simultaneous fit to Al ll A 1670 and Si ll A 1526 is used to 
derive the redshifts and b of the 4 components in the low- 
ionisation transitions of the sub-DLA at 2;abs=2.988 towards 
SDSS J0124+0400. These results are checked for consistency 
with many other Si ill lines (Si ll A 1193, Si ll A 1190, Si ll 
A 1260 and Si ii A 1526 - Si ii A 1808 is not covered by 
our data). The same redshifts and 6 are used to derive the 
column density of C ll A 1334 and upper limits from S ll A 
1253, S II A 1259 and Fe ill A 1122 which are all probably 
blended by lines from the Lyman-a forest. Similarly, an up- 
per limit is derived for Fe ll using the many lines falling in 
the forest (Fe ii A 1063, Fe ii A 1096, Fe ii A 1121, Fe ii A 
1125, Fe II A 1143 and Fe ll A 1144) as well as Fe ll A 1608. 
The fit is shown in Figure 5 and the matching parameters 
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Table 4. Parameters fit to the low-ionisation transitions of the Zabs=2.976 log 7Vhi=19-53±0.10 sub-DLA towards PSS J0121+0347. 



z b log N(Fe ll) 


log N(Si ii) 


log N{C II) 


log N{0 I) log N(A1 III) 


2.976613 9.90±0.94 13.18±0.30 


<13.15 


<13.72 


<13.98 11.93±0.32 



o.o 
1 .o 




CIV_1 548 



Ly forest 




Figure 4. Fit to the high-ionisation transitions of the 2:abs=2-976 
log7VHi=19.53±0.10 sub-DLA towards PSS J0121+0347 (see Ta- 
ble 5). 



Table 7. Parameters fit to the high-ionisation transitions of 
the Zabs=2.988 log ArHi=19.18±0.10 sub-DLA towards SDSS 
J0124-I-0044. 



z 


b 


log N(C iv) 


log N(Si IV) 


2.986426 


9.40±0.01 


13.60±0.02 


<13.00 


2.986677 


14.70±0.61 


14.19±0.18 


<13.59 


2.987084 


10.70±0.03 


13.44±0.02 


<13.34 


2.987509 


15.60±0.02 


13.31±0.08 


<13.41 


2.988232 


11.94±0.10 


13.10±1.04 


<13.34 


2.988896 


17.70±0.05 


13.01±0.07 


<13.47 


2.989483 


5.70±0.05 


12.33±1.41 


<12.97 




-100 -50 

Relative velocity (k-r/s) 

Figure 5. Fit to the low-ionisation transitions of the 2abs=2.988 
logAfHi=19.18±0.10 sub-DLA towards SDSS J0124-I-0044 (see 
Table 6). 



are presented in Table 6. The two AI ill lines are not covered 
by our spectrum for this system. 

The C IV doublet is nicely fitted with 7 components, the 
two first being the strongest ones. The Si iv doublet is 
blended but is nevertheless fitted using the same compo- 
nents in order to derive an upper limit on the abundance. 
Note however, that there is no physical reason suggesting 
that C IV and Si iv should come from the same region, i.e. 
should have the same velocity profile. It is just an empirical 
fact that in DLAs, this is often the case. The fit is shown 
in Figure 6 and the matching parameters are presented in 
Table 7. 



=3.840, z,bs=3.078, 



(iv) SDSS J0124-f-0044 
logArHi=20.21±0.10): 

This system has fairly high H l column density of 
logArHi=20.21±0.10. The Si ll A 1304 and Si ll A 1526 pro- 
files of the sub-DLA at Zabs=3.078 are well fitted by two 



components. The results of this simultaneous fit is applied 
on the blended Si ii A 1190, Si ii A 1193 and Si ii A 1260 for 
consistency check. The same redshifts and b are used to put 
limits on O I A 1302, C ll A 1334 and Fe ll A 1096, Fe ll A 
1125, Fe II A 1144 and Fe ll A 1608. O I and C ll are clearly 
saturated and lead to lower limits. A column density for Fe ll 
A is constrained from a combination of the red wing of Fe ll 
A 1608 and the blue wing of Fe ll A 1096 and is consistent 
with all the other Fe ll A lines available. The fit is shown 
in Figure 7 and the matching parameters are presented in 
Table 8. 

The Si IV doublet associated with this system is totally 
blended in the forest but the C iv seems to be detected as 
an extremely broad b>100 km s— 1 single component line. 
The strength of the two C iv lines of the doublet are con- 
sistent with their respective oscillator strengths, but they 
do not display the characteristic complex profiles of sub- 
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Table 6. Parameters fit to the low-ionisation transitions of the Zabs=2.988 log 7Vhi=19-18±0.10 sub-DLA towards SDSS J0124+0044. 



z 


b 


log N(Fe ii) 


log N(Si ii) 


log N{C II) 


log N(S II) 


log N(A1 II) 


log N(Fe III) 


2.986404 


2.90±0.01 


<12.62 


13.41±0.13 


<15.67 


<13.48 


11.68±0.16 


<13.39 


2.986680 


12.50±0.61 


<13.20 


13.98±0.11 


<14.62 


<13.59 


12.58±0.12 


<13.59 


2.988448 


12.80±0.03 


<13.00 


12.71±0.16 


<13.53 


<13.88 


11.82±0.13 


<13.51 


2.988998 


8.50±0.20 


<12.72 


12.78±0.16 


<14.07 


<13.63 


11.93±0.14 


<13.44 




Figure 6. fit to the high-ionisation transitions of the Zabs=2.988 
log7VHi=19.18±0.10 sub-DLA towards SDSS J0124+0044 (Ta- 
ble 7). 



DLA/DLAs. Therefore, care should be taken in interpret- 
ing these features as C iv. The column density we derive 
is N(C iv) = 13.92±0.11 but this system might as well be a 
large A^hi system with no C iv detected. This is illustrated 
in Figure 8. 

(v) PSS J0133+0400 (zcm=4.154, Zabs=3.139, 
logArHi = 19.01±0.10): 

This sub-DLA has a low A''hi column density with 
log A'hi = 19.01±0.10 at Zabs=3.139. Unfortunately, many of 
the expected strong metal lines for this system fall in spec- 
tral gaps (i.e. Fe ii A 1608, C iv A 1548, C iv A 1550) or 
are blended with H I lines from the Lyman-a forest (i.e. the 
other Fe ll A lines) . The two component fit was performed on 
Si II A 1260 with a consistency check on Si ll A 1808 to derive 
an upper limit on the total Sill column density. Attempts 
to increase the number of components does not improve the 
of the fit in this case. The same redshifts and b were ap- 
plied on the blended C ll A 1334 to derive an upper limit 
on the C abundance. The fit is shown in Figure 9 and the 
matching parameters are presented in Table 9. Al ill A 1854 
is strongly blended, but Al ill A 1862 is used to derive an 4(t 
upper limit: log N(A1 lll)<11.44. 

The position of the high-ionisation transitions are unfor- 
tunate in this system too. The Si iv A 1402 falls in a DLA 
along the same sight-of-line while the Si iv A 1393 is com- 
pletely blended with H I lines from the Lyman- a forest. The 
C IV doublet is situated in one of the spectral gap of our 
data. 
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Figure 7. Fit to the low-ionisation transitions of the 2abs=3.078 
logAfHi=20.21±0.10 sub-DLA towards SDSS J0124-I-0044 (see 
Table 8). 



(vi) PSS J0133-I-0400 (zom=4.154, Zabs=3.995, 
logAfHi = 19.94±0.15): 

This system has a fairly high A^'hi with 
logA'Hi=19.94±0.15 at 2abs=3.995. The redshifts and b of 
the 5 components were derived from partially blended lines 
of Si II A 1260 (telluric contamination) for two bluest com- 
ponents and Si ll A 1190 (Lyman-a forest contamination) 
for the three remaining components. The resulting fit was 
checked upon other Si ill lines (Si ll A 1193, Si ll A 1260 and 
Si II A 1808). As a result, we can obtain a good value of the 
total N(Si ll). Most of the Fe ll A lines in this system are 
blended with Lyman-a forest features (Fe ll A 1063, Fe ll A 
1096, Fe II A 1121 and Fe ll A 1125) or fall in a zero fiux 
gap (Fe II A 1143 and Fe ll A 1144) due to DLAs along the 
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Table 8. Parameters fit to the low-ionisation transitions of the Zabs=3.078 log 7Vhi=20.21±0.10 sub-DLA towards SDSS J0124+0044. 



z 


b 


log N(Fc II) 


log N(Si II) 


log N(C II) 


log N(0 I) 


3.077625 


6.20±3.20 


<13.95 


13.75±0.20 


>14.93 


>15.00 


3.077758 


2.10±0.40 


<13.65 


15.11±0.40 


>15.15 


>14.60 




Figure 8. High-ionisation ions of the absorber detected towards 
SDSS J0124+0044 with log ArHi=20.21±0.10 at 2abs=3.078. 



Table 9. Parameters fit to the low-ionisation transitions of 
the 2;abs=3.139 log AfHi=19.01±0.10 sub-DLA towards PSS 
J0133-I-0400. 



z 


b 


log N(Si II) 


log N(C II) 


3.138360 


31.60±0.20 


<12.87 


<13.25 


3.139605 


7.40±0.10 


<12.81 


<13.72 






1 862 
- 1 OO 



Figure 9. Fit to the low-ionisation transitions of the 2^1,5=3.139 
logAfHi=19.01±0.10 sub-DLA towards PSS J0133-I-0400 (see Ta- 
ble 9). The fit of C II A 1334 is used to derive an upper limit on 
the C II column density. 



Table 10. Parameters fit to the low-ionisation transitions 
of the 2abs=3.995 log AfHi=19.94±0.15 sub-DLA towards PSS 
J0133-I-0400. 



z 


b 


log N(Si II) 


3.994089 


3.10±0.16 


12.24±0.12 


3.994422 


8.50±0.11 


12.83±0.11 


3.995350 


5.90±0.21 


13.56±0.14 


3.995771 


20.00±0.15 


13.22±0.11 


3.996604 


4.70±0.11 


13.29±0.18 




Figure 10. Fit to the low-ionisation transitions of the 2abs=3.995 
logAfHi=19.94±0.15 sub-DLA towards PSS J0133+0400 (see Ta- 
ble 10). 



same line-of-sight. Only the Fe ll A 1608 line is covered by 
our spectrum but it is undetected. We derive an upper limit 
log N(Fe ll)<13.56 at 4a. Other ions are not detected and 
allow the derivation of upper limits: log N(Ni ll)<12.31, log 
N(A1 ii)<11.36 and log N(A1 iii)<11.69. The fit is shown 
in Figure 10 and the matching parameters are presented in 
Table 10. 

The 5 components for the high-ionisation transitions are 
derived from a simultaneous fit of the lines Si iv A 1393 and 
Si IV A 1402. Si iv A 1393 provides a column density deter- 
mination free from any contamination. The same redshifts 
and b are used to fit the C iv doublet of the system and 
derive an upper limit on this blended lines. The fit is shown 
in Figure 11 and the matching parameters are presented in 
Table 11. 
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Table 11. Parameters fit to the high-ionisation transitions 
of the Zabs=3.995 log 7Vhi=19.94±0.15 sub-DLA towards PSS 
J0133+0400. 



z 


b 


og N(C IV) 


log N(Si IV) 


3.994154 


8.90± 0.01 


<12.92 


12.99±0.25 


3.994482 


9.30± 0.86 


<12.94 


12.90±0.03 


3.994817 


9.30± 0.03 


<12.93 


12.47±0.49 


3.995580 


3.20± 0.92 


<11.95 


12.24±0.42 


3.996051 


15.60± 1.33 


<12.79 


12.89±0.01 




Relative velocity (krri/Ej) 

Figure 11. Fit to the high-ionisation transi- 
tions of the Zabs=3.995 log A''hi=19.94±0.15 sub-DLA towards 
PSS J0133-I-0400 (see Table 11). 



(vii) PSS J0133-h0400 (zcm=4.154, 2abs=3.999, 
logiVHi = 19.16±0.15): 

This system is at the low end of the column den- 
sity distribution for sub-DLAs with log A^hi=19.16±0.15 at 
Zabs=3.999. Most Si lll lines are covered in this absorber but 
many are blended with Lyman-a forest contamination (i.e. 
Si II A 1190, Si II A 1193) or telluric contamination (Si ll 
A 1260, Si II A 1526). Si ii A 1304 falls in a spectral gap. 
Therefore, the region around the expected position for Si ll 
A 1808 is used to derive an upper limit on the abundance 
of Si III: log N(Si ll)<14.12 at 4(t. Similarly, upper limits 
are derived for non-detection as follows: log N(A1 ll)<11.10 
from Al II A 1670, log N(A1 iii)<11.63 from Al iii A 1854 
and log N(Fe ll)<13.56 from Fe ll A 1608. 

For the high-ionisation transitions, a fit is performed on 
C IV A 1548 and C iv A 1550 simultaneously. The same 
redshifts and b are used for the Si iv A 1393/1402 doublet. 
All these lines appear broad and quite noisy and therefore 
lead to upper limits. The fit is shown in Figure 12 and the 
matching parameters are presented in Table 12. 

(viii) PSS J0133+0400 (2em=4.154, Zabs=4.021, 
logiVm = 19.09±0.15): 

This is also a system at the low end of the column density 
range of sub-DLAs with log A''hi=19.09±0.15 at 2abs=4.021. 
Again, many lines for this absorber are blended in the 
Lyman-Q forest. Upper limits from non-detections of Al ll A 



Table 12. Parameters fit to the high-ionisation transitions 
of the 2abs=3.999 log 7Vhi=19.16±0.15 sub-DLA towards PSS 
J0133-I-0400. 



z 


b 


log N(C IV) 


log N(Si IV) 


3.998718 


20.40±0.40 


<13.44 


<12.55 


3.999390 


34.60±0.30 


<13.47 


<13.15 




Relative v.loaity (km/,) 



Figure 12. Fit to the high-ionisation transi- 
tions of the 2abs=3.999 log 7Vhi=19.16±0.15 sub-DLA towards 
PSS J0133-I-0400 (see Table 12). 



1670, Si II A 1808 and Fe ll A 1608 are derived at 4cr level: log 
N(A1 ii)<11.20, log N(Si ii)<14.05 and log N(Fe ii)<13.56. 

A kinematically simple broad absorption is detected at 
the position of the Si iv doublet. The profiles are fitted using 
both lines simultaneously to derive an upper limit. The C iv 
doublet however is not detected in this system. An upper 
limit on the abundance is derived to be log N(C iv)<12.31. 
The fit is shown in Figure 13 and the matching parameters 
are presented in Table 13. 

(ix) BRI J0137 4224 (zcm=3.970, Zabs=3.101, 
logiVm = 19.81±0.10): 

The lower redshift sub-DLAs of this quasar has many 
metal lines falling in the Lyman-a forest. With a column 
density log A'^hi=19.81±0.10 at 0abs=3.1Ol, it is situated in 
the mid-range of the sub-DLA definition. Nevertheless, some 
strong lines are undoubtfuUy clear from any contamination. 
Si II A 1526 and Al ll A 1670 in particular are used to de- 
rive the redshifts and b of the 7 components used to fit the 
metal-line profile. The resulting z and b are used to derive 
the abundance of Fe ll from a fit to the Fe ll A 1608 line. The 



Table 13. Parameters fit to the high-ionisation transitions 
of the 2abs=4.021 logArHi=19.09±0.15 sub-DLA towards PSS 
J0133-I-0400. 



z b log N(Si IV) 

4.021325 30.3±1.00 <12.89 
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Relative veleeity (krr,/;.) 



Figure 13. Fit to the high-ionisation transi- 
tions of the ^abs=4.021 log7VHi=19.09±0.15 sub-DLA towards 
PSS J0133+0400 (see Table 13). 




0.0 I , , t , , , , t , , , , 1 , , , , 

-50 O 50 100 

Relative veiocity (km/5) 

Figure 14. Fit to the low-ionisation transitions of the 2:abs=3.101 
logArHi=19.81±0.10 sub-DLA towards BRI J0137-4224 (see Ta- 
ble 14). 

remaining Fe ll A lines are heavily blended with Lyman-a 
forest contamination (i.e. Fe ll A 1096, Fe ll A 1144). The 
fit is shown in Figure 14 and the matching parameters are 
presented in Table 14. 

Concerning the high-ionisation transitions, the Si IV dou- 
blet is falling in the Lyman-a forest, therefore totally 
blended and not useful for abundance determination. Simi- 
larly, the C IV doublet is blended with unidentified lines as 
well as with apparent emission lines which probably are the 
products of bad cosmic rays clipping (no telluric lines were 
identified in this region) and have been removed to perform 
the fit. The C iv A 1550 is fitted with 8 components and the 
resulting profile is made consistent with C iv A 1548. The 
fit is shown in Figure 15 and the matching parameters are 
presented in Table 15. 

(x) BRI J0137 4224 (2om=3.970, Zabs=3.665, 
log7VHi = 19.11±0.10): 



Table 15. Parameters fit to the high-ionisation transitions 
of the 2abB=3.101 log AfHi=19.81±0.10 sub-DLA towards BRI 
J0137-4224. 



z 


b 


log N(C IV) 


3.100324 


3.50±0.16 


12.74±0.18 


3.100453 


3.20±0.07 


12.43±0.10 


3.100640 


4.50±0.28 


12.55±0.14 


3.100995 


3.10±0.24 


12.38±0.10 


3.101195 


20.40±0.12 


12.34±0.34 


3.101350 


3.20±0.23 


12.57±0.10 


3.101517 


3.20±0.10 


12.79±0.25 


3.101620 


9.90±0.20 


12.83±0.13 




Figure 15. Fit to the high-ionisation transi- 
tions of the 2abs=3.101 log AfHi=19.81±0.10 sub-DLA towards 
BRI J0137-4224 (see Table 15). 

This sub-DLA has a fairly low column density of 
log Ai'Hi=19.11±0.10 at 2;abs=3.665. Nevertheless, it shows 
strong, kinematically simple, metal lines which are well fit- 
ted by a single component. Three lines (O I A 1302, Si ll 
A 1260 and C ll A 1334) clear from any contamination are 
fitted simultaneously to derive the redshift, b and appropri- 
ate column densities. The resulting fit is overplotted on the 
many blended Si ill lines (Si ll A 1304, Si ll A 1190, Si ll A 
1190, Si II A 1526 and Si ii A 1808) and is found to be con- 
sistent in all cases. Unluckily, many of the Fe ll A lines for 
system fall in regions contaminated by Lyman-a forest lines 
(Fe II A 1063, Fe ii A 1096, Fe ii A 1121, Fe ii A 1125, Fe ii A 
1143 and Fe ll A 1144). Only Fe ll A 1608 lies in a region free 
from any contamination but the S/N of the spectrum is low. 
This leads to an upper limit log N(Fe ll)<13.59 at 4(t. Al ll 
A 1670 and Al lll A 1854/Al lll A 1862 are covered by our 
spectrum but no lines are detected. We derive upper limits 
on the column densities of these ions: log N(A1 ll)<11.06 
and log N(A1 lll)<11.96. The fit is shown in Figure 16 and 
the matching parameters are presented in Table 16. 

The high-ionisation transitions of this system present a 
very particular signature. Both lines of the Si iv doublet 
(1393 and 1402) are clearly detected and very well fitted 
by a simple two component profile. But rather surprisingly, 
both C IV lines (1548 and 1550) of the C iv doublet are 
covered by our data but seem undetected. In fact, the weak 
component falling at the expected redshift in C iv A 1550 
happens to be a telluric line, although the C iv A 1548 region 
is free from telluric contamination. At any rate, this line is 
within the noise of the spectrum and could not be assumed 
to be a real line for sure. From the S/N in the region where 
the doublet is expected to fall, we deduce an upper limit on 
the column density log N(C iv)<12.11 at 4a. Interestingly in 
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Table 14. Parameters fit to the low-ionisation transitions of the 2aba=3.101 log 7Vhi=19-81±0.10 sub-DLA towards BRI J0137-4224. 



z 


b 


log N{Fe II) 


log N{Si II) 


log N(A1 III) 


3.100382 


2.90±2.64 


12.65±0.39 


12.92±0.33 


11.73± 0.14 


3.100559 


4.90±0.11 


12.95±0.01 


13.40±0.20 


11.88± 0.89 


3.100821 


7.40±0.91 


12.70±0.01 


13.64±0.29 


12.20± 0.72 


3.101004 


3.30±0.12 


12.80±0.12 


12.73±0.22 


11.40± 0.54 


3.101679 


8.00±0.77 


13.20±0.01 


13.45±0.25 


12.33± 0.49 


3.101915 


3.20±0.22 


12.70±0.01 


13.25±0.10 


12.06± 0.10 


3.102078 


4.10±0.40 


12.10±0.01 


12.27±0.20 


11.48± 0.20 




Figure 16. Fit to the low-ionisation transitions of the 2:abs=3.665 
logAfHi=19.11±0.10 sub-DLA towards BRI J0137-4224 (see Ta- 
ble 16). The region 25-25 km s~^ is contaminated by an uniden- 
tified blend. 




Relative velocity (km/s) 



Figure 17. Fit to the high-ionisation transi- 
tions of the ^abs~3. 665 logArHi=19.11±0.10 sub-DLA towards 
BRI J0137-4224 (see Table 17). 



Table 17. Parameters fit to the high-ionisation transitions 
of the Zabs=3.665 log AfHi=19.11±0.10 sub-DLA towards BRI 
J0137-4224. 



z 


b 


log N(Si IV) 


3.665009 


6.60±5.20 


12.43±0.22 


3.665482 


7.30±3.80 


12.95±0.10 



Dessauges-Zavadsky et al. (2003) , we report the detection in 
a sub-DLA of C iv lines with no Si iv associated. But to our 
knowledge, this is the first time that a large A^hi absorber is 
reported to have Si iv but no C iv detected. This is could be 
a signature of Si iv associated with the A^hi gas whilst C iv 
will be part of an external shell surrounding this region. The 
fit is shown in Figure 17 and the matching parameters are 
presented in Table 17. 

(xi) BR J2215 1611 (^0^=3. 990, Zabs=3.656, 
logA'Hi = 19.01±0.15): 

The two sub-DLAs detected towards BR J2215— 1611 are 
essentially blended together and only come apart at the 
Lyman-7 level. The first system at 2abs=3.656 is also the 
weakest H l column density with log AfHi=19.01±0.15. Metal 
lines are seldom detected partly because of contamination 



with the Lyman-a forest lines. Few lines (Si ll A 1304, Si ll 
A 1526, Si II A 1808 and Fe ll A 1608) faU in regions free 
from any contamination and yet remain undetected. We 
derive upper limits on the column densities of the corre- 
sponding ions: log N(Si ll)<14.41, log N(Fe ll)<13.51 and 
log N(S ll)<13.68. The intermediate ionisation transitions 
Al III A 1854 and Al lll A 1862 are falling in the region pol- 
luted by telluric lines which prevent us from any detection 
or even determination of an upper limit. 

The high-ionisation transitions, C iv A 1548 and C iv A 
1550, are detected in this system. To the immediate blue 
wavelength of the C iv A 1548 line lies a telluric line. The 
profile is well fitted by one component whose characteristics 
are derived from a simultaneous fit of the two members of 
the doublet, but only an upper limit is derived to be conser- 
vative. The Si IV in this sub-DLA is not detected although 
both lines are covered by our data. This leads to an up- 
per limit on the column density: log N(Si iv)<11.69. The 
fit is shown in Figure 18 and the matching parameters are 
presented in Table 18. 

(xu) BR J2215 1611 (zcm=3.990, Zabs=3.662, 
logAfHi=20.05±0.15): 

Contrary to the previous system, this sub-DLA has 
a fairly high A'hi column density log A^hi=20.05±0.15 at 
-2abs=3.662. It is characterised by strong metal lines show- 
ing a complex profile over --^100 km s~^. The 5 components 
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Table 16. Parameters fit to the low-ionisation transitions of the 2aba=3.665 log A'Hi=19.11it0.10 sub-DLA towards BRI J0137— 4224. 



z 6 log N(Si II) log N(C II) log N(0 i) 

3.665111 7.02±0.50 12.40±0.13 13.14±0.13 13.38±0.13 




0,5 I , I 1 , 1 , I 1 , , I 1 , . 1 I 1 , 

—SO O SO 1 OO 

R^lotive velocity (km/s) 

Figure 18. Fit to the high-ionisation transi- 
tions of the -Sabs— 3. 656 logAfHi=19.01±0.15 sub-DLA towards 
BR J2215-1611 (see Table 18). 

are simultaneously fitted to O I A 1302, the un-saturated 
components of C ii A 1334, Al ii A 1670 and Si ii A 1526. 
The first component is narrower in Si ll A 1526 than in the 
other metal lines including Si ll A 1304. We note however 
that O I A 1302 is controversial and that it is a borderline 
case between a lower limit and a value. 

The resulting fit is applied to the other Si ll lines: Si ll A 
1304 is clearly weU fitted by this profile. Si ll A 1808 is not 
detected in agreement with our fit, while Si ll A 1190, Si ll 
A 1193 and Si ll A 1260 which are situated in the Lyman-a 
forest are clearly blended by other intervening lines. Nev- 
ertheless, the fit is found to be consistent with the upper 
hmit they provide. All Fe ll A lines but Fe ll A 1125 are 
contaminated, but we use the non detection of Fe ll A 1125 
to derive log N(Fe ll)< 13.60. The fit is shown in Figure 19 
and the matching parameters are presented in Table 19. The 
intermediate ionisation transitions Al ill A 1854 and Al ill A 
1862 are falling in the region polluted by telluric lines which 
prevent us from any detection or even determination of an 
upper limit. 

The high-ionisation transitions in this system are well fit- 
ted with 3 components. The redshifts and b of these were 
determined by a simultaneous fit of the C iv A 1548 and 
Si IV A 1402 lines. The fit is then superimposed on C iv 



Table 20. Parameters fit to the high-ionisation transitions 
of the ^abs=3.662, log AfHi=20.05±0.15 sub-DLA towards BR 
J2215-1611. 



z 


b 


log N(C IV) 


log N(Si IV) 


3.661409 


7.20± 5.00 


12.79±0.07 


12.67±0.08 


3.661881 


12.40± 3.10 


13.57±0.25 


13.39±0.04 


3.662255 


26.50± 8.80 


12.92±0.08 


12.53±0.07 



A 1550 and Si iv A 1393. Clearly the Si iv A 1393 line is 
blended although no telluric lines are seen in this region. 
This fit leads accurate column density determinations. The 
fit is shown in Figure 20 and the matching parameters are 
presented in Table 20. 

(xiii) BR J2216-6714 (zcm=4.469, Zabs=3.368, 
logArHi=19.80±0.10): 

This sub- 
DLA has a mid-range column density log A^hi = 19.80±0.10 
at Zabs=3.368. Nevertheless, the low-ionisation transitions in 
this system are barely detected, partly because of Lyman-a 
contamination (the metal lines with Arest <1522 A in this 
sub-DLA will fall in the forest) . This applies to many of the 
Si III lines (Si ii A 1304, Si ii A 1260, Si ii A 1193, Si ii A 
1190 and Si ii A 1526). The non-detection of Si ii A 1808 lead 
to the determination of an upper limit: log N(Si ll)<13.88 
at 4a. Similarly, most of the Fe II lines are in the forest, 
but the region at the expected position of Fe ll A 1608 has 
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Table 19. Parameters fit to the low-ionisation transitions of the 2aba=3.662, log A'^Hl=20.05it0.15 sub-DLA towards BR, J2215— 1611. 



b log N(Si ii) log N(C ii) log N(0 i) log N(A1 ii) 



3.661105 
3.661597 
3.661754 
3.661885 
3.662193 



3.30±0.10 
6.80±0.10 
6.80±0.20 
5.10±0.10 
3.50±0.10 



12.95±0.28 
13.45±0.02 
12.60±0.03 
13.38±0.01 
13.11±0.02 



13.72±0.48 
>13.95 
>13.35 

14.81±0.20 
>14.16 



14.03±0.01 
14.61±0.20 
13.61±0.02 
14.25±0.01 
14.58±0.02 



11.70±0.16 
11.92±0.27 
11.47±0.10 
11.95±0.67 
11.55±0.53 




Relative velocity (km/.) 



Figure 20. Fit to the high-ionisation transi- 
tions of the Zabs=3.662, log AfHi=20.05±0.15 sub-DLA towards 
BR J2215-1611 (see Table 20). 

a fairly high signal-to-noise leading to a constraining upper 
limit on the Fe ll column density: log N(Fe ll)<13.26. Al ill 
A 1854 and 1862 are also undetected: log N(A1 lll)<11.51. 
On the contrary, a two component feature at the expected 
position of Al ll A 1670 is clearly detected. We compute an 
abundance determination for this element. The fit is shown 
in Figure 21 and the matching parameters are presented in 
Table 21. 

In addition, the high-ionisation transitions are also un- 
doubtfuUy detected for this system. The C iv doublet is 
clearly detected and leads an accurate column density de- 
termination using a simultaneous fit of both members of the 
doublet to determine the redshifts and fo of 6 components. 
The Si IV doublet appears to be blended with Lyman-a for- 
est interlopers, although Si iv A 1393 is most probably free 
form any contamination. The parameters derived from C iv 
are used to determine the upper limits on Si IV A 1393, but 
the first and fourth components are not detected in Si IV. 
The fit is shown in Figure 22 and the matching parameters 
are presented in Table 22. 



3 RESULTS 

3.1 Individual Abundances 

The resulting total column densities for each elements de- 
tected in the 13 sub-DLAs studied here are summarised in 
Table 23 together with the error estimates. 

One concern is that sub-DLAs may be partially ionised 



Table 21. Parameters fit to the low-ionisation transitions 
of the Zabs=3.368, log AfHi=19.80±0.10 sub-DLA towards BR 
J2216-6714. 



z 


b 


log N(A1 II) 


3.369173 


15.20±2.54 


11.88±0.26 


3.370053 


17.40±3.04 


11.85±0.27 



Figure 21. Fit to the low-ionisation transi- 
tions of the 2j,bs=3.368, log ArHi=19.80±0.10 sub-DLA towards 
BR J2216-6714 (see Table 21). 

in H, in which case the observed Fe ll to H I ratio would 
not be a good measure of the total Fe abundance for in- 
stance. Indeed, given that Fe ll has an ionisation potential, 
IP(FeII) = 16.18eV, larger than hydrogen, IP(HI) = 13.59eV, 
some of the observed Fe ll could reside in the ionised gas of 
the sub-DLA with lower A^hi. This leads to a slight overesti- 
mate of the true total metallicity of the systems. To investi- 
gate the ionisation corrections, we compute photo-ionisation 
models based on the CLOUDY software package (version 
94.00, Ferland 1997), assuming ionisation equilibrium and 
a solar abundance pattern. Based on the modelling of 13 
systems, Dessauges-Zavadsky et al. (2003) have shown that 
differences between a Haardt-Madau or a stellar-like spec- 
trum are negligible. Here, we chose to perform the modelling 
using a stellar-like incident ionisation spectrum (D'Odorico 

Table 22. Parameters fit to the high-ionisation transitions 
of the Zabs=3.368, logArHi=19.80±0.10 sub-DLA towards BR 
J2216-6714. 



z 


b 


log N(C IV) 


log N(Si IV) 


3.368701 


20.10± 4.07 


12.95±0.11 




3.369287 


15.50± 2.96 


13.02±0.12 


<12.52 


3.369919 


12.60± 5.42 


13.01±0.09 


<12.34 


3.370177 


6.10± 1.99 


12.58±0.22 


* 


3.370583 


74.00± 1.59 


13.49±0.08 


<13.09 


3.370990 


8.30± 2.01 


13.03±0.08 


<12.48 



*: These components are undetected in Si IV. 
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Table 23. Total column densities for each absorber. 



Quasar ^^bs log N(Fe ii) log N(Si ii) log N(C ii) log N(0 i) log N(A1 ii) log N(A1 in) log N{C iv) log N(Si iv) 



PSS J0118+0320t 


4.128 


14.16±0.11 


14.84±0.13 


>17.30 


>16.40 






14.30±0.12 


13.78±0.12 


PSS J0121+0347 


2.976 


13.18±0.30 


<13.15 


<13.72 


<13.98 


11.93±0.32 


<11.77* 


<13.79 


13.33±0.10 


SDSS J0124+0044t 


2.988 


<13.55+ 


14.12±0.12 


<15.72 




12.76±0.13 




14.43±0.30 


<14.20 




3.078 


<14.13 


15.13±0.39 


>15.35 


>15.15 










PSS J0133+0400 


3.139 




<13.14 


<13.85 






<11.44* 








3.995 


<13.56'® 


13.91±0.13 






<11.36* 


<11.69* 


<13.51 


13.48±0.28 




3.999 


<13.56* 


<14.12* 






<11.10* 


<11.63* 


<13.76 


<13.25 




4.021 


<13.56* 


<14.05* 






<11.20* 




<12.31* 


<12.89 


BRI J0137-4224 


3.101 


13.67±0.11 


14.11±0.25 






12.83±0.42 




13.52±0.17 






3.665 


<13.59* 


12.40±0.13 


13.14±0.13 


13.38±0.13 


<11.06* 


<11.96* 


<12.11* 


13.06±0.13 


BR J2215-161lt 


3.656 


<13.51' 


<14.41* 










<13.22 


<11.69* 




3.662 


<13.60* 


13.89±0.06 


>14.98 


15.05±0.02 


12.46±0.45 




13.71±0.21 


13.51±0.05 


BR J2216-6714 


3.368 


<13.26* 


<13.88* 






12.17±0.26 


<11.51* 


13.88±0.10 


<13.32 



*: 4 (T upper limits corresponding to non-detections 
+ : log N(Fe lll)<14.09 in this system 

®: log N(Ni ll)<12.31 in this system from a non-detection 

t: PSS J0118+0320, 2a,bs=4.128 has log N(S ii)=14.26±0.11, SDSS J0124-f0044, Zabs=2.988 has log N(S ii)<14.27, BR 
J2215-1611, Zabs=3.656 has log N(S ll)<13.68 from a non-detection. 




O 50 1 OO 1 50 200 250 300 



Relative velocity (km/s) 

Figure 22. Fit to the high-ionisation transi- 
tions of the Zabs=3.368, log AfHi=19.80±0.10 sub-DLA towards 
BR J2216-6714 (see Table 22). 



& Petitjean 2001). We thus obtain the theoretical column 
density predictions for any ionisation state of all observed 
ions as a function of the ionisation parameter U. When pos- 
sible, we then use observed ratios of the same element in 
different states (i.e. Fe ll/Fe ill, Al ll/Al ill. Si lll/Si iv or 
C ll/C iv) and compare them with predictions from the 
model to constrain the U parameter. When only limits on 
the U parameter were available, we have established a limit 
on the ionisation correction. Moreover, on two occasions, 
such ratios are not available and the observed column densi- 
ties are directly compared with the predicted ones to deduce 
the U parameter. 

Once the ionisation parameter U is known, we can de- 
duce by how much the observed metallicity, Fellobs/ Hlobs, 
deviates from the total metallicity, Fctot /Htot of the ab- 
sorber. We compute the correction to bring to the observed 
Fe II values, the so-called Sioni as follows: 




18.5 19 19.5 20 20.5 



loq N(HI) 

Figure 23. Ionisation correction, Sioni, as a function of logA^jji 
from the present study and from the literature (Prochaska et 
al. 1999; Dessauges-Zavadsky et al. 2003; Peroux et al. 2006a; 
Prochaska et al. 2006). There seems to be trend of smaller ionisa- 
tion corrections towards higher A^hi, as expected. The solid line 
is the linear least squares regression with slope a = 0.13. Note 
that in most cases, the correction, Sioni, is smaller than 0.2 dex. 

_ Fctot Fellobs 

Oioni — „ „ , (Ij 

tlTOT tllobs 

A negative (positive) Sioni corresponds to an overes- 
timate (underestimate) of the total metallicity. Figure 23 
shows the correction for 26 sub-DLAs as function of A'hi 
column densities from the present study and from the litera- 
ture (Prochaska et al. 1999; Dessauges-Zavadsky et al. 2003; 
Peroux et al. 2006a; Prochaska et al. 2006) . Recent findings 
by Fox et al. (2007) on the ionised fraction in O and C of a 
sample of DLAs are also in line with the present results. The 
solid line is a fit to the measures (as opposed to lower limits) 



Global Metallicity Evolution 13 



• • 

# 


_ 


• t 




• • 






• 







-1 

Metollicity 



Figure 24. lonisation correction, Sioni, as a function of observed 
metallicity, [Fe/H] (or [Si/H] when Fe measurements are not avail- 
able). There is a hint of decreasing Sioni with metallicity. 




Absorption Redshift 

Figure 25. lonisation correction, &ioni, as a function of absorp- 
tion redshift. Given than the incident UV background flux evolves 
with redshift, a correlation of the lonisation corrections with z^^,^ 
is expected. No such correlation is observed. 



with slope a=0.13. As in Dessauges-Zavadsky et al. (2003) 
and more recently, Prochaska et al. (2006), we find that the 
lonisation correction for sub-DLAs are small and within the 
error estimates in most cases (sse also Erni et al. 2006 for a 
borderline case) . In the present study, only 4 sub-DLAs out 
of the 13 studied require a correction |(5ioni| > 0.2 dex. In all 
these 4 cases, the correction is \5ioni\ < 0.35 dex. The Sioni 
values for each sub-DLAs are listed in Table 24. There is a 
hint of decreasing Sioni with metallicity (Figure 24). Know- 
ing that sub-DLAs are more metal-rich than classical DLAs 
(Kulkarni et al. 2007), possibly due to an observational bias, 
and that the lonisation correction increases with decreasing 
A'hi (see Figure 23); it is no surprise to observe that the loni- 
sation correction gets larger as the metallicity gets higher. It 




Redshift 

Figure 26. Redshift evolution of [Fe/H] metallicity in DLAs 
and sub-DLAs. Filled circles are sub-DLAs studied by our group 
(Dessauges-Zavadsky et al. 2003 and the present study). Sub- 
DLAs arc found to evolve more with time than DLAs, and to be 
more metal-rich at the lowest redshifts. The lines show the least 
square regressions fitted to the data (excluding limits) for DLAs 
(dashed line) and sub-DLAs (solid line). 



is interesting to note however that the fractional correction 
remains roughly constant. However, no clear correlation of 
Sioni with absorption redshift is detected (Figure 25), con- 
trary to what one would expect from the redshift evolution 
of the incident UV background photons. 

The total absolute abundances were calculated with re- 
spect to solar using the following formula: 



[X/H] = log[N{X)/N{H)]i 



log[N{X)/N{H)]Q (2) 



where log[N{X)/N{H)]Q is the solar abundance and 
is taken from Asplund et al. (2005) adopting the mean of 
photospheric and meteoritic values for Fe, Si, C, O, S and 
Al. These values are recalled on the top line of Table 24. 



3.2 Global Metallicity 

Figure 26 shows the resulting redshift evolution of the metal- 
licity of DLAs and sub-DLAs including our new sample of 
13 high-redshift sub-DLAs. The data plotted are based on 
Peroux et al. (2003b) compilation with the following studies 
added: Dessauges-Zavadsky, et al. (2004), Rao et al. (2005), 
Peroux et al. (2006a), Peroux et al. (2006b) and Meiring et 
al. (2007). We note that thanks to these recent studies the 
metallicity of sub-DLAs is now better constrained. These 
results confirm findings from Peroux et al. (2003b) that the 
average values of the metallicities in the difi'erent redshift 
bins do show more evolution with redshift in the case of 
sub-DLA than for classical DLAs. Fitting least square re- 
gressions to the data points without including limits lead 
to a slope of a = —0.18 for DLAs and a — —0.42 for sub- 
DLAs. These regression fits are plotted as dashed and solid 
lines (respectively) on Figure 26. Moreover, based on Zn 
measurements, sub-DLAs appear more metal-rich at lower 
redshifts, specially at z < 2 (Peroux et al. 2006a; Peroux 
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Table 24. Abundances with respect to solar, [X/H], using the standard definition: [X/H] = loglN(X)/N{H)]nLA - log[N{X)/N{H)]Q . 
The error bars on [X/H] include both the errors in log N{X) and logAfjji- Sioni is the correction to the observed metaUicity required to 
take into account the ionised part gas (see equation 1). 
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et al. 2006b, Meiring et al. 2007). Although our data add 
no new Zn measurements because of their high-redshifts, it 
has been previously shown (Peroux et al. 2003b) that the 
strong evolution observed in [Fe/H] for sub-DLAs is not due 
to differential depletion given that such strong evolution is 
also observed in [Zn/H] metaUicity of sub-DLAs. 

Rather, the difference in evolution might be explained 
by the fact that sub-DLAs are less prone to the biasing effect 
of dust and thus represent a better tool to detect the most 
metal-rich galaxies seen in absorption. Indeed, it has been 
suggested that the dust contained in the absorbers might 
introduce biases into current surveys in the sense that the 
dustier systems would dim the hght from the background 
quasar and therefore not be taken into account in current 
magnitude-limited samples. In order to overcome such limi- 
tations, searches for DLAs towards radio-selected quasars 
have been undertook and their results suggest that dust 
obscuration might be modest (Ellison et al. 2001; Jorgen- 
son et al. 2006). However, the samples are still small to de- 
rive firm conclusions and some of the radio-selected quasars 
are actually "dark" at optical wavelengths and can thus po- 
tentially be obscured quasars. In fact, if dust extinction in 
quasar absorbers is a strong function of N(Zn ll) column 
density as found in interstellar clouds (Vladilo et al. 2004), 
the obscuration will start acting at a lower [Zn/H] ratio for 
DLAs than for sub-DLAs (Vladilo & Peroux, 2005). In other 
words, we might be missing more systems in the DLA range 
than in the sub-DLA range and the metal-rich sub-DLAs 
recently found would be the "tip of the iceberg" population 
which has remained so far un-noticed. A way to find metal- 
rich absorbers minimising dust effects is therefore to search 
for quasar absorbers with relatively lower A^'hi column den- 
sity such as the sub-DLAs. Another possibility to explain 
the metaUicity difference between DLAs and sub-DLAs has 
recently been suggested by Khare et al. (2006). Based on 
the observed mass-metallicity relationship for galaxies (e.g. 
Tremonti et al. 2004; Savaglio et al. 2005), the authors argue 
that sub-DLAs may arise in massive galaxies and DLAs in 
less massive galaxies. 




19 20 21 22 

log HI 

Figure 27. [Fe/H] metaUicity of both DLAs and sub-DLAs as 
a function of log A^hi . The dotted vertical line corresponds to 
the DLA definition: log Niii> 20.3. Different redshift ranges are 
symbol-j-colour coded and show once more the fastest redshift- 
evolution of sub-DLAs with respect to DLAs. 

In Figure 27, the metaUicity of both DLAs and sub- 
DLAs is shown as a function of log A^'hi. The spread in metal- 
licity is larger in the sub-DLAs range showing once more a 
strong evolution with redshift, while classical DLAs appear 
more homogeneous. 



4 DISCUSSION 

4.1 The Missing Metals Problem 

A direct consequence of the star formation history is the 
production of heavy elements, known as metals. However, 
at high-redshift, our knowledge of the cosmic metal bud- 
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get is still highly incomplete. In fact, the amount of metals 
observed in high-redshift galaxies (e.g. Lyman-Break Galax- 
ies, DLAs) and the intergalactic medium, was believed to be 
roughly a factor of 5 below the expected amount of metals 
produced as a result of the cosmic star formation history. 
This is dubbed as the "missing metals problem" (Pettini 
1999 and Pettini 2004). This census has recently been revis- 
ited in a series of papers, which explore alternative popula- 
tions which might also contain part of the missing metals. 
Bouche, Lehnert & Peroux (2005) have investigated the role 
of sub-mm galaxies while Bouche, Lehnert & Peroux (2006) 
have considered the Lyman-Break Galaxies. But a substan- 
tial fraction of the missing metals may also be hidden in 
a very hot, coUisionally ionised gas (Ferrara et al. 2005). 
Based on simple order-of-magnitude calculations, Bouche et 
al. (2007) discuss the possibility that the remaining missing 
metals could have been ejected from small galaxies via galac- 
tic outflows into the intergalactic medium in a hot phase 
which is difficult to detect using observed properties of local 
galaxies. However, even when taking into account the most 
recent observations, it appears that 10 to 40% of metals are 
still missing. 

Although we can measure the metallicities of DLAs up 
to high-redshifts, the above studies find that their ad-hoc 
metallicities is too low for them to be a major contribu- 
tor to the metal budget. However, Dave & Oppenheimer 
(2007) used cosmological hydrodynamic simulations based 
on GADGET-2 to compute the redshift evolution and con- 
tribution of DLAs to the metal census. They find that their 
simulations reproduce well the mild redshift evolution ob- 
served in DLAs but overproduce the metallicity at all red- 
shifts. They suggest that highly enriched sub-DLAs might 
accommodate for the discrepancy. In the following, we use 
the most recent observational data on sub-DLAs to calculate 
their contribution to the metal budget. 

4.2 Contribution of Sub-DLAs to the Global 
Metallicity 

Recently, there have been several attempts to estimate the 
amount of metals contained in sub-DLAs (Prochaska et al. 
2006), including some results based on dust-free Zn metallic- 
ities measures (Kulkarni et al. 2006). Indeed, Pettini (2006) 
has suggested that some of the missing metals might be in 
sub-DLAs if those were to have a significant ionised frac- 
tion. But the amount of metals in the ionised gas of sub- 
DLAs cannot be directly probed by observations. In fact, 
it is interesting to note that Lebouteiller et al. (2006) have 
recently shown that the ionised gas probed in emission is 
more metal-rich that the neutral gas probed in absorption 
for most elements in a local HII region. While the reasons 
for this is not yet clear, they also note that Fe ll is a better 
tracer of neutral gas than Si ill. 

For high-redshift quasar absorbers, the only mean to 
quantify the ionised fraction is to use photo-ionisation mod- 
els. The ionised fraction is defined as: 

X = /Htot 

Prochaska et al. (2006) estimate a:=0.90 for one sub- 
DLA and compute Q.z{sub — DLAs) using this value. Here, 
we use the models presented in Section 3.1, to reproduce the 
ionisation state of 26 sub-DLAs. The x values for these are 
shown in Figure 28 as a function of A^hi column densities. 




loq N(HI) 

Figure 28. Ionised fraction, x = H'^ / HxOT^ of a sample of 26 
sub-DLAs as a function of A^jji column densities. The up and 
down arrows denote lower and upper limits, respectively. Note 
that a number of x values are found at any given A'^jjl column 
density (see for example AfHi~19.7). 
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Figure 29. Evolution of the ionised fraction, x, with the metal- 
licity, [Fe/H] (or [Si/H] when Fe measurements are not available). 
No obvious trend is observed. 



A range of x values are found at any given A^hi column 
densities. Figure 29 shows the evolution of the ionisation 
fraction with metallicity for the 26 sub-DLAs studied. No 
obvious trend is observed. 

At z ~2.5, the amount of HI gas in sub-DLAs is mea- 
sured to be r^Hi^O.lS X 10"^ (Peroux et al. 2005). At the 
same redshifts, the metallicity of these systems is measured 
from the undepleted Zn element is Z/Zq=— 0.70 (Kulkarni 
et al. 2007). The mean of a;-values (as opposed to limits) at 
z ^^2.5 is < a; >= 0.68. Therefore, assuming all the gas is 
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Figure 30. Evolution with cosmic time of the ionised fraction, 
X, in sub-DLAs. 



photo-ionised, the observed comoving mass density of metals 
in sub-DLAs is: 



xl0"°'^°x0.18xl0"^x 



overdensities at different redshifts. In fact, systems which 
are more affected by higher ultra-violet flux probably have 
a lower observed A^hi and are therefore not classified as sub- 
DLAs. Indeed, the gas cooling time is expected to be ~ 10^ 
yrs, i.e. much longer than the time it takes for H I to turn 
into stars (~ 10^ yrs) or to be heated by the meta-galactic 
flux. This also goes in line with the fact that only small 
amounts of molecular hydrogen H2 are found in quasar ab- 
sorbers (Srianand et al. 2005, Zwaan & Prochaska 2006). 



5 SUMMARY 

We have presented abundance and ionisation fraction deter- 
minations of a sample of 13 2: > 3 sub-DLAs. In summary: 

• our new high-resolution observations more than double 
the metallicity information for sub-DLAs previously avail- 
able at 2: > 3. 

• results from photo-ionisation modellings of a sample of 
26 sub-DLAs are put together for the first time. 

• the ionisation correction to the observed metallicity, 
Sioni is a function of the observed A^'hi column density and 
smaller than 0.2 dex in most cases. 

• the metallicity of sub-DLA evolves faster with cosmic 
times than the one of classical DLAs and is higher especially 

-al(3) < 2. This could be due to sub-DLAs being less affected 



nz{sub - DLAs) = 2.57 x 10" 



(4) 



where ^©=0.0126 by mass and f2(Zo) — ilbaryons x 
Zq = 5.5 X 10"''. 

The total amount of metals expected at this redshift 
is Q{Z) = 4.5 X 10~^. Therefore, the sub-DLAs are con- 
tributing at most ~ 6% of the total metal in the Universe at 
z ~ 2.5 for a photo-ionised gas. According to recent results, 
if the gas was collisionally ionised, the contribution to the 
total amount of metals would be similar (Fox et al.2007). 
So, even if sub-DLAs are mostly ionised gas, they do not 
close the metal budget. In fact, given that 10 to 40% of the 
metals are still missing, the mean metallicity of sub-DLAs 
after ionisation correction should be —0.46 < Z/Zq < +0.14 
for these systems to close the missing metals problem. Such 
metal-rich galaxies would be easily observed in emission. 
This is not the case of quasar absorbers which probably 
trace fainter objects. Note that the higher the ionised frac- 
tion of sub-DLAs, X, the higher is their contribution to the 
cosmic metallicity. It will thus be very important in the fu- 
ture to measure the ionised fraction of more sub-DLAs, to 
better constrain < x > at z ~2.5. 



4.3 Redshift Evolution of the Ionised Fraction 

Figure 30 shows the evolution with cosmic time of the 
ionised fraction in sub-DLAs. Clearly, one would expect the 
ionised fraction of sub-DLA to increase with cosmic times as 
the Universe becomes ionised and the meta-galactic flux in- 
creases up to z=1.5. On the contrary, no clear trend seems 
to appear with redshift. This result supports the picture 
where sub-DLAs are a phase of galaxy evolution observed 
at various redshifts rather than a tracer of the continuous 
formation of galaxies and/or that we are probing different 



by the biasing effect of dust. 

• the ionisation fraction of a sample of 26 sub-DLAs al- 
lows to reliably calculate the total contribution of these sys- 
tems to the so-called "missing metals problem". Sub-DLAs 
contribute no more than 6% of the total amount of expected 
metals at z'^2.5. 
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